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We study static annihilation on complex networks, in which pairs of connected particles annihilate 
at a constant rate during time. Through a mean-field formalism, we compute the temporal evolution 
of the distribution of surviving sites with an arbitrary number of connections. This general formal- 
ism, which is exact for disordered networks, is applied to Kronecker, Erdos-Renyi (i.e. Poisson) 
and scale-free networks. We compare our theoretical results with extensive numerical simulations 
obtaining excellent agreement. Although the mean- field approach applies in an exact way neither 
to ordered lattices nor to small-world networks, it qualitatively describes the annihilation dynamics 
in such structures. Our results indicate that the higher the connectivity of a given network element, 
the faster it annihilates. This fact has dramatic consequences in scale-free networks, for which, once 
the "hubs" have been annihilated, the network disintegrates and only isolated sites are left. 
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I. INTRODUCTION 

The reaction dynamics of particles that undergo pair- 
wise mutual annihilation is of interest to various fields 
of physics, including diffusion controlled reactions^, Q , 
and exciton annihilation in molecular crystals [3j. In 
these examples the particles are free to move or diffuse 
throughout the system, and when two complementary 
particles collide, they annihilate. This kind of dynam- 
ics can be modeled by adopting a network approach in 
which the particles, rather than moving across the sys- 
tem, are placed at the vertices of a given network. The 
edges of the network would represent direct interactions 
between the particles at a given time, or interactions that 
might occur during time as the particles diffuse through- 
out the system. This approach has proven to be use- 
ful in understanding the collective dynamical behavior in 
many-particle systems 0, 0- I n some other cases, the 
annihilating particles can be accurately represented as 
being fixed on a lattice. As an example, consider the py- 
cnonuclear reactions taking place in the interior of dense 
stars [a, til ISl- In this case, the atom nuclei undergoing 
mutual annihilation become frozen into a regular crys- 
talline lattice structure due to the high densities of the 
cores. Other examples are the irreversible sequential ad- 
sorption and "car parking" analogues, where the system 
reaches a jammed state in which the lattice or space can- 
not be filled completely 

In the present paper we investigate the static pairwise 
annihilation problem on complex networks, gene ralizing 
previous results obtained for ordered lattices jlOl Y\A\ . As 
such, the annihilation problem on complex networks is 
interesting in its own right since it provides a dynami- 
cal probe of the topological properties of the structure 
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on which it occurs. These properties will be different 
than, and complementary to, those associated with other 
physical processes in complex networks, such as random 
wal ks I f 21 . Ising spin dynamics 1% target annihilation 
|l4l Il5j, majority voter models and neural network dy- 
namics 0, ^3 1 epidemic spreading 0, 0] , trapping 
[HE1, etc. 

It is precisely these complementary features which lead 
us to consider static annihilation as a probe of different 
complex networks, such as ordered lattices, Erdos-Renyi 
random networks (also called Poisson networks), scale- 
free networks and small- world networks. Small- world 
networks are particularly interesting, since they can be 
considered as a class of networks lying at some interme- 
diate point between ordered lattices and fully random 
networks. They possess a large degree of local cluster- 
ing, but a relatively small minimal path length connect- 
ing nodes throughout the system. On the other hand, 
scale-free networks are characterized by the presence of 
"hubs" or elements with an extremely high connectivity. 
These hubs hold the network together and produce a very 
small minimal path length. Complex networks are cur- 
rently under intense investigation as useful structures on 
which to model a variety of network related phenomena, 
including electrical power grids, the Internet, human so- 
cial relationships, proteomic and genetic interactions, to 
mention just a few examples [53, HjJ 0> HE US • We are 
interested here in exploring the effects that arise on the 
static annihilation process from the different features as- 
sociated with the different network topologies mentioned 
above. 

Several exact results have previously been obtained 
for the static annihilation problem on ordered lattices. 
Kenkre and Van Horn, e.g., considered a ID ordered ring 
of N sites, in which each site is initially occupied by a 
particle, as in Fig. [IJi. At the onset of the annihilation 
process, neighboring particles on the ring begin to un- 
dergo mutual annihilation with rate constant a. As pairs 
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FIG. 1: Schematic diagram illustrating the annihilation pro- 
cess in a ID ring, (a) At time t = there are N particles 
placed along the ring (black circles), (b) As time goes by, 
pairs of neighboring particles annihilate at a constant rate a. 
At some time t > some pairs have been annihilated (open 
circles). Note that all the neighbors of the particle signaled 
with an arrow have been annihilated. Since this particle has 
no neighbors left to annihilate with, it will survive indefinitely. 

of particles disappear from the ring, the possibility arises 
for particles at certain sites to become isolated. Such 
particles end up surviving indefinitely, with no neighbors 
left on either side to annihilate with (see Fig. |Ip) . For 
this model the authors of Ref. 0>E1 were able to calcu- 
late exactly the average fraction fit) of sites on the ring 
that remain populated at time t. For an infinite ID ring 
they found 

f(t) = exp (2e' at ~ 2) , (1) 
from which the exact stationary survival fraction 

= lim f(t) = e" 2 (2) 

t — >oo 

follows. A dynamical mean-field theory for higher dimen- 
sional ordered structures was later given by Kenkre 
In that work, the long-time survival fraction on a square 
lattice in D dimensions was found to decrease as the di- 
mensionality D (and the number of potential annihilation 
partners) increases, specifically, via 

f oc = (2D-l) D ^ 1 - D \ (3) 

Of course, on a translationally invariant structure the 
connectivity of any site on the lattice is the same as that 
of any other. Each site, moreover, has a priori the same 
probability of ultimately surviving the annihilation on- 
slaught. It is interesting, therefore, to see what effects 
are generated in the annihilation dynamics by dispersion 
in the local site coordination number. In a sense, Pois- 
son, scale-free and small-world networks are ideal struc- 
tures on which to explore these questions. Such struc- 
tures allow, through a simple variation of the network 
parameters, significant changes in the average coordina- 
tion number as well as a means to systematically change 
the associated dispersion in local coordination. This is es- 
pecially true for scale-free networks, where the dispersion 
in the local coordination number can become infinite. 




(a) (b) 

FIG. 2: (a) A fc-site (at the center) with k connecting sites. 
k = 8 in this diagram, (b) One of the connecting sites of the 
fc-site is an m-site with m connections (m = 4 in the diagram). 
The probability that an arbitrary connection (dashed line) of 
the fc-site goes to an m-site is mf(m, t)/ X^m=o m f{m,t). 



In the present paper we generalize the results of |10llll| 
to a variety of complex networks that allow us to explore 
these connectivity-related issues. First, in Section 1TT1 we 
present a mean-field theory of the annihilation process 
on disordered networks with arbitrary connectivity dis- 
tributions. We then apply the general formalism to three 
specific cases: (a) random networks in which every ele- 
ment has exactly K connections (Kronecker networks); 
(b) Erdos-Renyi (i.e. Poisson) networks with average 
connectivity K; and (c) scale-free networks. It is worth 
mentioning that the mean-field approach requires the sta- 
tistical independence of the network elements. Thus, 
our formalism applies neither to ordered lattices nor to 
small-world networks, where loops of connected sites are 
commonly found. Such loops generate non-trivial corre- 
lations between the elements which are not taken into 
account in the mean-field approach. However, case (a) 
above can be considered as an approximation to ordered 
lattices with coordination number K. The approxima- 
tion becomes exact in the limit K — > oo. Additionally, 
as we show through numerical simulation in Section IlIII 
the mean-field approach qualitatively describes the anni- 
hilation dynamics in ordered and small-world structures. 
We end this work in Section IIVI with a summary and 
discussion of our results. 



II. MEAN-FIELD THEORY 

In this section we develop a simple mean-field theory of 
particle annihilation in networks with an arbitrary con- 
nectivity distribution P{k). This theory is exact as long 
as the connecting sites of each element in the network are 
randomly chosen from anywhere in the system. In such a 
case, the correlations that arise due to the occurrence of 
loops of connected sites can be neglected. The above is 
not true in general either for small-world networks or for 
ordered lattices. However, the mean-field approach pre- 
sented here will allow us to gain insight into the effects 
that the overall connectivity of the system has on the 
annihilation dynamics. As we will see, since the infinite 
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ID ordered lattice with nearest-neighbor connections is 
loopless, our predictions recover the exact solution for 
this structure. 

We are interested in evaluating the fraction of surviv- 
ing sites f(t) at time t in a network in which dimers (i.e. 
pairs of connected sites) have been removed randomly. 
We assume that the network is initially full, and that the 
probability that a site has fc connections is given by P(k). 

We denote by f(k, t) the fraction of surviving sites that 
are connected to k surviving sites at time t. We will refer 
to these as fc-sites. It is clear that f(k,t = 0) = P(k). 
Next, we determine the time evolution of these densities 
by considering the three different ways in which /(fc,t) 
can change: 



The above equations can be treated by generating func- 
tions as follows. Let F(z,t) be the generating function 
defined as 



F(z,i) = ]Tz fc /(M). 



(•5) 



Note that /(i) = F(l,t) = Efcto /(M) is J ust the overall 
fraction of surviving elements at time t. From Eq. Q it 
follows that F(z, t) satisfies the equation 



dFM =[{l- Z ) 9 {t)- Z ] dF ^ 



dt 



dz 



(6) 



1. As we remove dimers, a fc-site can become a (fc — 1)- 
site if one of its connecting sites is removed. To de- 
termine the rate at which this happens we focus on 
one of the k connections of the fc-site (dashed line 
in Fig. |2t>) ■ The probability that this connection 
goes to an m-site is proportional to mf(m,t) |3^ . 
We have to normalize this probability by dividing 
between $Zm=i m /( TO ; ^) since the fractions /(m, t) 
change in time and they do not add up to 1 (ex- 
cept at t = 0). Therefore, the chosen connecting 
site in Fig. |3) will be an m-site with probability 

mf(m,t)/ J2m=i m f( m 't)- The m-s he will be re- 
moved if the dimer formed by this site and any of 
the other (m — 1) sites to which it is connected (not 
considering the fc-site we started from) is removed. 
Adding over all m and considering the fc connec- 
tions of the fc-site, the total rate at which fc-sites 
become (fc — l)-sites is 



fc/(fc,f) 



Em=O m /(«M) 



2. A (fc + l)-site can become a fc-site. For the same 
reasons as before, this happens with probability 



(fc + l)/(fc + M) 



E^o m ( m - l)f(m,t) 



3. Finally, fc-sites may also be removed if any of the fc 
different dimers to which they belong is chosen for 
removal, which occurs with probability fc/(fc, t). 



where 



9(t) 



m(m — l)f(m,t) 

m=0 

oo 

J2 mf(m,t) 

m=0 



( d 2 F(z,t) 

f dF(z,t) 
\ dz 



(7) 



z=l 



Equation © can be solved by the method of character- 
istics, which gives: 



F(z,t) 
where 



l + (z-l)e 



-G(t) 



and 



G(t) = t+ g{r)dT, 
Jo 



Q(z) = 5> fc P(fc). 

fc=0 



(8) 



(9) 



(10) 



Note that Q(z) is the generating function of the initial 
distribution P(k). It is also important to note that g(t) 
must be determined through the consistency equation 
0. 

It is clear that at the end of the annihilation process, 
only isolated elements with fc = will survive. Therefore, 
the fractions f(k,t) should fulfill the condition 



Adding the three quantities above (with the correct 
sign) we find that the densities f(k,t) satisfy the follow- 
ing set of evolution equations: 



df(k,t) 
dt 



= [(fc + l)/(fc + l,t)-fc/(fc,i)] x 



J2 m(m ~ l)/(m,t) 

rn— 

oo 

E mf(m,t) 

m=0 



-kf(k,t). (4) 



lim/(M) = fora11 k > 1 - 

t — >oo 



The fraction of surviving elements at the end of the 
annihilation process is then given by 



f 00 = lim f(0,t)= lim F(l,t). (11) 

t — >oo t — >oo 



In what follows we apply the general formalism pre- 
sented in this section to obtain for networks with 
different topologies. 
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A. Kronecker Delta Distribution 

We now apply the preceding formalism to the specific 
case in which every site of the initial network has exactly 
K connections. Therefore, the initial distribution P(k) is 
given by a Kronecker delta function: P(k) = 8k,K- This 
case can be considered as an approximation to an ordered 
lattice with coordination number K. Such an approxi- 
mation becomes exact for K — » oo. However this case 
is particularly interesting when K = 2 because, given 
that in the "thermodynamic" limit there are no loops of 
connected sites in the network, the system corresponds 
exactly to the infinite ID lattice with nearest-neighbor 
connections, for which the exact solution is known (see 
Eqs. OJ and © and Refs. QUEl). For this choice 01 
the initial distribution, Eq. (|10|l becomes 



Q(z) = 5> fc P(fc) 



k=0 



From this result we find that the generating function for 
the distribution of fc-sites, Eq. (JSJ, is given by 



r r* 1* 

F{z, t)= 1 + [z - l)e- G W - / e- G(T) dT 

Jo 

and the consistency equation (Q) becomes 
g(t) = -{K - l)j t In (l - £ e' G ^dr 



(12) 



(13) 



To solve the preceding integral equation we first inte- 
grate it from to t, which gives 



/ g(r)dT = -(K-l)ln\l- f 
Jo l Jo 



or, equivalently, 

cxp - / g{r)dT 
. Jo 



= 1- 



t \ K-l 

e- G ^dr 



Multiplying both sides of the above equation by e 'we 
obtain 



e -G(t) =e -tl 1 



e~ G ^dr 



(14) 



where we have used the definition of G(t) given in Eq. 
The following change of variable will be very useful to 
solve the integral equations: 



u{t) 



e- G ^dr. 



(15) 



With this change of variable, Eq. I|14|l can be written as 
du(t) 



dt. 



[«(*)] 



K-l 



(16) 



The above equation can be solved easily. However, 
care must be taken for the case K = 2, which must be 
evaluated separately. Assuming K > 2 first, the solution 
of Eq. I|16f) with the initial condition u(0) = 1 is 



i(t)=[l + (K-2){l-e- t )] 



_ t sil/(2-K) 



(17) 



Although the complete hierarchy of distributions 
f(k, t) can be obtained from the preceding results, at this 
point we are interested in the final fraction of surviving 
sites /oo. This can be readily computed by noting from 
Eqs. and JTSJ) that the overall fraction f(t) = F(l, t) 
of surviving sites at time t is 



/(*) = [«(*)] 



K 



[l + (A--2)(l-e- 



t ,iK/(2-K) 



Taking the limit t — * oo in the above equation we obtain 



foc = (K— l) K/{ 



2-K) 



(K > 2). 



(18) 



Since a square lattice in D dimensions has coordination 
number K = 2D, the above result coincides with Eq. J2J, 
which was first reported in Ref. |l 1| . 

As mentioned before, the case K = 2 has to be evalu- 
ated separately. For K = 2, the solution of Eq. I|16(l with 
the initial condition u(0) — 1 is 



u(t) = exp (e * - l) 



(19) 



Therefore, the overall fraction f(t) of surviving sites in 
this case becomes 

/(*) = [u{t)f = exp (2e-* - 2) , 

which, in the limit t — > oo, gives 

/co = e- 2 (K = 2). (20) 

This is the exact result for the ID lattice with nearest- 
neighbor connections. (See Eq. (J2J and Refs. [Tol ITl| 

Fig. 01 shows /oo as a function of K . The solid curve is 
the theoretical result as predicted in Eqs. (|TH|l and lt20ll . 
whereas the circles correspond to the computer simula- 
tion data. Note from Eq. Ijl8|l that /oo ~ 1/K for large 
values of K, and therefore 

Jim /oo = 0. 

So, the more connected the initial network, the smaller 
the fraction of surviving sites at the end of the annihila- 
tion process. 

For the sake of completeness, we compute the explicit 
form of the distributions f(k, t) for the case K > 2. To 
this end, we rewrite Eqs. I|15(l and (|16|l as 

e-GW = e - t [u(t)] K - 1 , 



f e- G{T) dr = l-u(t). 
Jo 
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FIG. 3: Fraction of surviving elements foo as a function of the 
average network connectivity K for networks with Kronecker 
distribution. The circles are the plot of the simulation data 
whereas the solid curve is the graph of the theoretical predic- 
tion given in Eqs. 1181 and 12011 . The dashed line indicates the 
value e -2 corresponding to the ID ring with nearest-neighbor 
interactions. To generate the simulation data we used net- 
works with TV = 10 3 sites. Each point is the average over 100 
realizations. 



Substituting the above expressions into Eq. (|12J) we ob- 
tain for the generating function 

F(z,t)= [(z-l)e- t [u(t)] K - 1 + u(t)] K , 

which, after expanding in powers of z using the binomial 
theorem, transforms into 



A (K\ (uffl-e-Xt)]*- 1 ) 
fe W {et[u{t)?-K-lf 



K 



From this result it follows immediately that the dis- 
tribution f(k,t) of surviving sites with k connections at 
time t is given by 

/(M) = (f) ( ° (t) - C ~' Ktr T. PD 

W (e<HI)p-*' - 1)' 

where for u(t) we have to use Eq. H17fl if K > 2, or 
Eq. (fH3|) if K = 2. In any case, from Eq. (|2"TJ) we can 
see that lim^oo f(k, t) = for all fc > 1. Namely, at the 
end of the annihilation process, only isolated sites with 
k = survive, as expected. However, for k > 1 the rate 
at which f(k,t) — > as t — > oo depends on the value of k. 
To see this, consider the distributions f(k,t) and f(k',t) 
with k > k'. Then, from Eq. (|2"TJ> we get 



hm 



lim 



t^oo f(k>, t) t^oo ( e t[ M (t)]2-K _ lf- k ' 



0. 



This result shows that more connected sites annihilate 
faster than less connected ones, as expected. 



B. Poisson Distribution 

Next we treat the case in which the initial distribu- 
tion of connections is Poisson with mean K: P(k) = 
e~ K K k /k\. This case corresponds to the Erdos-Renyi 
topology, for which Q{z) becomes 



Q { z)=Y J Z k P{k)=e- K Y.z k ^=z- K{1 

k=0 fc=0 

Using this result in Eq. © we obtain 



(22) 



F(z, t) = exp 



-if (1 - z)e 



and the consistency equation (JJJ in this case is 



g(t) = Kcxp i-t - J g{r) drX = Ke 



G(t) 



(23) 



(24) 



The above integral equation can be easily solved for g(t), 
obtaining 



g(t) 



Ke 



l + if(l-e-*)' 



(25) 



From Eq. I|24|) we have e G ^ — g(t)/K, from which it 
follows that 

P -G(t) _ e _ 



I + K{1 - e-ty 
1 



e~ G ^dT = j-lnll + Kil-e-*)}. 

Substituting the above results into Eq. l(2"3*|) we obtain 

e -g(t) 



F(z,t) 



1 + K(1- e-*) 
e -g{t) 



1 + K(l - e~*) ^ fc! 

v ; fe=0 



E 



(26) 
(27) 



where we have used the fact that g(t) — Ke~ G ^\ It 
follows immediately from Eq. (|27|1 that the distribution 
f(k,t) of fc-sites at time t is 



/(M) = 



,-g{t) 



W)] k 



1 + K(1 



k\ 



(28) 



An interesting aspect of this result is the fact that the 
distribution of fc-sites continues to be Poissonian (though 
not normalized) throughout the evolution of the system. 
Note also that, since lim^oo g(t) = (see Eq. P5fr). 
it follows from Eq. that lixa t ^oo f(k,t) — for all 
k > 1, as expected. Again, the rate at which this happens 
depends on the value of the connectivity k. Indeed, con- 
sidering the distributions f(k, t) and f(k', t), with k > k', 
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FIG. 4: Fraction of surviving elements /oo as a function of 
the average network connectivity K for Erdos-Renyi networks 
with Poisson distribution. The circles are the plot of the simu- 
lation data whereas the solid curve is the graph of the theoret- 
ical prediction given in Eq. 1291 . Each point of the simulation 
data is the average over 100 realizations using networks with 
N = 10 3 sites. 



from Eq. I|28|) and the fact that g(t) — > as t — > oo, we 
obtain 



,. /(M) 
-7(M 



lim [g(t)} 

t — >oc 



k-k' 



0. 



Therefore, the higher the connectivity of a given site, the 
faster it annihilates. 

The overall fraction f(t) of surviving sites at time t, 
regardless of their connectivity, is obtained from Eq. I|26|) 
with z = 1, which gives 



f(t)=F(l,t) = 



1 



l + K(l-e- f )' 



Finally, the fraction /oo of surviving sites at the end of 
the annihilation process can be obtained from the above 
equation (or from Eq. I|28ll with k — 0) by taking the 
limit t — > oo. We obtain 



1 



(29) 



1 + K 

Fig. 0] shows the graph of as a function of K. The 
solid curve is the plot of the theoretical prediction given 
in Eq. 1)29(1 and the circles correspond to the simulation 
data. Note that, once again, 
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as for the Kronecker-delta case. 
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C. Scale-free distribution 

In this section we consider scale- free networks for which 
the initial distribution of connections acquires the form 



P(k) 



1 



C(7) 



where C(t) — SfcLi ^ s the R-iemann zeta function. 
In this case, the generating function 1(1 0|l of the initial 
distribution becomes 



sv " k=l 



k=l 
1 



Li 7 (z), 



(30) 



where Li 7 (z) = X^feLi z k /k~ f is the polylogarithm func- 
tion of order 7. This function has the following proper- 
ties, which will be of use for further analysis: 



Li 7 (0) = 0; 

Li 7 (l) = C(t); 
dLi-tjz) 1 

— = -Ll 7 _i(2) 

oz z 



(31a) 
(31b) 

(31c) 



<9 2 Li 7 (z) 1 



dz 2 



— [Li 7 _ 2 (z)-LVi(z)]. (31d) 



The generating function F(z, t) for this case is 

F(z, t) = ^Li 7 (l + (z- l)e- G ^ - J e-^dr^j , 

(32) 

and, using properties l(31b ) and H31t 0. the consistency 
equation J7J becomes 

g (t) = e' GW Li 7 _ 2 ( M (t))-Li 7 _!( M (t)) 



u(t) 



where u(t) is defined as in Eq. (|15fl . Note from Eqs. 
and 1(15(1 that the overall fraction f(t) of surviving parti- 
cles at time t is given by 



f(t)=F(l,t) 



C(7) 



Li 7 (u(t)) ■ 



(34) 



Therefore, to compute f(t) we first have to know u(t). In 
order to do so, we rewrite the consistency equation 133|) 

as 

ff (t) = lpn(«(t))-ln(Li 7 _ 1 («(t)))], 

where we have used the fact that du(t)/dt = —e~ G ( l \ 
The last equation can be integrated from to t taking 
into account that u(0) — 1, which leads to 

' g{t)dt = In ( ^— ) + In (« 7 - 1)) - 
\Li7-l («(*))/ 

Adding t on both sides of the previous equation and ex- 
ponentiating we get 

-G(t) _ Li 7 _i(u(f)) e~* 



«(*) ' C(7 " 1) 



or equivalently, 



u(t) dtt(f) 



LL_i(u(t)) dt C(7-l)' 



(35) 
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We have been so far unable to find an exact analytic so- 
lution to the above nonlinear differential equation. How- 
ever, we can solve it numerically to obtain the asymp- 
totic value «oo = lim^oo u(t) with any desired accuracy. 
Once Moo has been determined, it follows from Eq. I|34|) 
that the final fraction /oo of surviving particles at the 
end of the annihilation process will be given by 



foo 



C(7) 



LLy (««,). 



(36) 



In order to determine Woo, we first integrate Eq. 1|35|) 
from to t, which gives 



du' = 



(t) 1%-iK) C(7-!) 



Taking the limit t — > oo in the previous equation we 
obtain 



«oo Li 7 -i( u ') 



-du' 



1 



C(7 - 1) 



(37) 



An important conclusion can immediately be drawn from 
this equation. Note that, since C(7 — 1) = oo for 1 < 7 < 
2, the above equation in this range becomes 



-du' = for 1 < 7 < 2, 



which has the solution u^ = 1. From Eq. i|36|) and prop- 
erty l|31bi. we obtain 

f OQ = l for 1< 7 < 2. 

In other words, essentially all the particles of a scale- 
free network with infinite average connectivity survive the 
annihilation process. This result is easy to interpret. 
What holds a scale-free network together are the "hubs" , 
i.e., the sites with an extremely high number of connec- 
tions But the hubs represent a negligible fraction 
of the total number of sites in the network. However, 
a randomly chosen site will be connected to a hub with 
certainty. Therefore, in this range of values of 7, the 
hubs are annihilated from the very beginning of the pro- 
cess. But once the hubs have been annihilated, the net- 
work disintegrates and only isolated sites are left. Hence, 
/oo = 1 in the thermodynamic limit. 

To obtain /oo for 7 > 2 it is convenient to define the 
function / 7 (u) as 



I 7 (u) 



1 u' 



-du'. 



' Li 7 i{u'Y 
With this definition, Eq. 137f) can be written as 

1 



(38) 



(39) 



This is a transcendental equation for . Note that J~ (u) 
is a continuous monotonically decreasing function of u. 



lyfa) 0,4 




FIG. 5: Features in scale-free networks, (a) The solid curve 
is the graph of J-y(tt) (defined in Eq. 1381 ') as a function of 
u for the case 7 = 2.5. The dashed line indicates the value 
1/£(1.5) ~ 0.382793, whose intersection with I~,(u) in this 
case gives w 0.431. (b) as a function of 7. This 
curve was obtained by solving numerically the transcendental 
equation 13911 . 



For example, the solid line in Fig. [SJl is the graph of 
Ij(u) for 7 = 2.5, whereas the dashed line indicates the 
value of 1/C(7 — 1), which for 7 — 2.5 is 1/£(1.5) s» 
0.382793. Then, Uqo is the value of u at which the solid 
curve intersects the dashed one, which for the case shown 
in Fig. is Moo ~ 0.431. Since I 7 (u) is a "well behaved" 
function for every 7, we can find the roots of Eq. (|39() with 
any desired accuracy. Fig. [3Jd shows Uoo as a function of 
7. Using these results in Eq. (|36|) . we obtain the data 
reported in Fig.|SJi, in which /oo is plotted as a function of 
7. The solid curve depicts the theoretical result and the 
circles are data from the simulation of the system in the 
computer. The small difference between the theoretical 
prediction and the simulation at 7 = 2 is due to finite 
size effects in the computer simulation. 

Finally, Fig. |SJd shows the graph of /oo as a function 
of the average network connectivity K, which for scale- 
free networks is given by K = £(7 — 1)/C(7)- Note that 
/oo increases with K and asymptotically approaches its 
maximum value /oo = 1 as K — > 00. As explained before, 
due to the presence of "hubs" in scale-free networks, we 
have 



Urn /oo 

K— *oo 



contrary to what happens in Kronecker and Poisson net- 
works, for which lim^-^oo /oo = 0. 
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FIG. 6: (a) Surviving fraction f x in scale-free networks as a 
function of 7. The solid curve corresponds to the theoretical 
result as predicted by Eqs. 1361 and 13711 . whereas the circles 
are the data obtained from computer simulation. The simu- 
lation was carried out for networks with N = 10 4 sites. Each 
point is the average over 100 realizations, (b) Same data as 
before but plotted against the average network connectivity 
JT = C(7-1)/C(7)- 



III. EXTENSION TO SMALL- WORLD 
NETWORKS 

Small-world networks pose a special challenge for ana- 
lytical treatment since we have to deal with both ordered 
short-range and random long-range interactions simulta- 
neously. Although some analytic results have been ob- 
tained for percolation and ferromagnetic-like systems in 
small- world networks p?l l28l 12^ . 130 ] , to the best of our 
knowledge there does not yet exist a general formalism to 
study the dynamical properties of this kind of networks 
given a specific interaction rule between the elements. In 
particular, for the annihilation problem we are consider- 
ing, the ordered short-range connections produce "loops" 
of connected sites, which in turn give rise to nontrivial 
correlations between the elements. Such correlations are 
not taken into account in the mean-field theory we have 
presented. Therefore, we do not expect our formalism to 
be exactly applicable to small-world networks. 

However, since small- world networks can be considered 
as something between ordered lattices and fully random 
networks, we do expect the static annihilation process 
in SW structures to be similar to that in Kronecker or 
Poisson networks. In this section we present results ob- 
tained from the numerical simulation of the system on 
small-world networks. Such networks are constructed by 
starting with an ordered ring of N sites and connectiv- 
ity K in which each site is connected to its K/2 near- 
est neighbors on the left and K/2 nearest neighbors on 



FIG. 7: (a) Surviving fraction f x in small-world networks as 
a function of K for three values of the rewiring probability. 
These three values correspond to an ordered ring (p = 0), a 
small- world network with 20% of shortcuts (p = 0.2), and a 
fully random network (p — 1). (b) Fraction of surviving el- 
ements foa as a function of p for three values of the average 
network connectivity K. In both figures the data were com- 
puted numerically through computer simulation using net- 
works with N = 10 3 sites. Each point is the average over 
1000 realizations. 



the right. Disorder is then introduced into this ordered 
structure by disconnecting with probability p each bond 
in the lattice emerging from a given site and reconnecting 
it to another randomly chosen site to which it was not 
originally connected (provided that the rewiring does not 
result in two or more disconnected subnetworks). 

Thus, the reconnection or rewiring probability p is a 
measure of the disorder of the system, and allows the 
system to be tuned from a completely ordered (p = 0) 
to a completely disordered (p = 1) structure. It is also a 
direct measure of the number of "short-cuts" that exist 
in the network. Note that the rewiring algorithm leaves 
the average connectivity K associated with the original 
ordered structure unchanged, but changes the dispersion 
in the local connectivity. Indeed, after rewiring any given 
site may be connected to a greater or fewer number of 
sites than before rewiring. Therefore, the average con- 
nectivity and the magnitude of fluctuations in the local 
connectivity are independently adjustable through the 
parameters K and p. 

Fig-EK shows foo as a function of K for three different 
values of the rewiring probability: p = (ordered ring), 
p = 0.2 (small-world network) and p = 1 (fully random 
network). Note that behaves as ~ 1/K for large 
values of K regardless of the value of the rewiring prob- 
ability p. Therefore, small-world networks behave like 
Kronecker and Poisson networks. On the other hand, 
Fig. [7jp shows foo as a function of p for three values of 
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the average network connectivity K. Interestingly, the 
numerical results of this figure indicate that the frac- 
tion /oo of surviving particles increases with an increas- 
ing number of small-world shortcuts. One might think 
that such shortcuts would enhance particle annihilation 
between sites very far away, and thereby lead to a de- 
crease in the survival fraction. However, in order for 
a small-world network to achieve the short path length 
with which it is associated, a significant number of sites 
have to be more connected than others in the system. 
Such sites act as gateways, allowing access across the 
system from any direction. Since the average network 
connectivity is preserved during rewiring, the existence of 
highly connected elements implies the existence of many 
elements with a connectivity smaller than average. Con- 
sequently, there will exist situations in which sites with 
high connectivity are connected to many sites of below 
average connectivity. As our mean- field formalism shows, 
a highly-connected site is very likely to annihilate quickly 
with one of its many partners, leaving the remaining sites 
to which it was connected more isolated than they would 
have been if they had been part of an ordered network. 
Hence, it appears that dispersion in the local connec- 
tivity leads to enhanced survivability in the annihilation 
problem. 

IV. SUMMARY AND DISCUSSION 

In this work he have analyzed the static annihila- 
tion process on complex networks. In this problem, 
pairs of connected particles annihilate at a constant rate. 
Through a mean-field approach we were able to com- 
pute exact general expressions for the fraction of sur- 
viving particles with k connections at time t for ran- 
dom networks with arbitrary connectivity distribution 
P(k). Our results indicate that Kronecker and Pois- 
son networks exhibit a similar behavior in that the fi- 
nal fraction of surviving particles /oo asymptotically be- 
haves as 1/K for large values of K, the average con- 
nectivity of the initial network. Numerical simulations 
show that this is also true for ordered lattices and small- 
world networks. So, although the mean-field approach is 
not exactly applicable to ordered latices and small-world 
networks, it predicts qualitatively the correct asymptotic 
limit /oo ~ 1/K for such systems. 



The situation is radically different in scale-free net- 
works, characterized by the connectivity distribution 
P(k) ~ fc~ 7 . For this kind of network, our theory pre- 
dicts that /oo = 1 for 1 < 7 < 2 and that Iim 7 _ +O o /oo 
1 . Thus, in the range 1 < 7 < 2 in which the aver- 
age network connectivity K diverges, all the particles (or 
almost all in finite systems) survive the annihilation pro- 
cess. The above result is a consequence of the fact that 
highly connected sites annihilate faster than poorly con- 
nected ones. While this happens in every network, it has 
dramatic consequences in scale-free structures in which, 
once the hubs arc annihilated, the network breaks apart 
and only isolated sites are left. This kind of behavior is 
reminiscent of the breakdown of scale- free networks when 
the highly connected sites are deliberately attacked [T3 |. 
Here we have shown that the annihilation of the hubs 
does not only break the network apart, but that it pul- 
verizes the network. 

Diffusion annihilation dynamics of the kind A + A — > 
has recently been studied on complex networks using an 
analogous formalism as the one presented in this work 
|3lj |. In diffusion annihilation, the particles, rather than 
being fixed to the vertices, are free to move (diffuse) 
throughout the network. Two particles annihilate if they 
collide at some vertex of the network for the first time. 
The authors of Ref. [jyj show that the fraction of sur- 
viving particles decreases in time as t~ x for Poisson-likc 
networks and as t~@ for scale-free networks, where f3 de- 
pends on the scale- free exponent. Obviously, contrary to 
what happens in the static case considered in this work, 
there are no surviving particles in diffusion annihilation. 
Thus, our results describe situations that are, in some 
sense, complementary to those analyzed by Ref. j2j| ■ 
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